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Using 4 flavors of staggered fermions at infinite gauge coupling, we compare var- 
ious analytic results for the hadron spectrum with exact Monte Carlo simulations. 
Agreement with Ref. is very good, at the level of a few percent. 

Our results give credence to a discrepancy between the baryon mass and the 
critical chemical potential, for which baryons fill the lattice at zero temperature 
and infinite gauge coupling. Independent determinations of the latter set it at 
about 30% less than the baryon mass. One possible explanation is that the nuclear 
attraction becomes strong at infinite gauge coupling. 
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1 Introduction 



The strong coupling limit of lattice QCD can provide valuable insight, as Wilson 
showed by proving the area-law behaviour for the Wilson loop in that regime [Tj. 
Analytic integration over the gauge link variables becomes feasible in this limit. 
Within some reasonable approximation (typically mean field), one can then calcu- 
late analytically the hadron spectrum. The result turns out to be not too differ- 
ent from the continuum QCD spectrum. This has been the subject of numerous 
works 121 El in E] • To the best of our knowledge however, these analytic results 
involving different approximation schemes have never been compared with the re- 
sults of numerical simulations in the infinite coupling limit (/5 = 6/(7^ = 0). Here, 
we perform such a comparison for the case of 4 flavors of staggered fermions with 
SU(3) gauge group. 

In addition, we are motivated by a puzzle in the phase diagram of QCD as a 
function of temperature T and chemical potential fi. This phase diagram also can 
be determined analytically when (3 = 0, and has been investigated by use of increas- 
ingly sophisticated approximations over the years jSj-CI]- At zero temperature, a 
first-order transition is predicted to take place, where the baryon density jumps 
from zero to saturation (1 baryon per lattice site). The corresponding critical 
quark chemical potential for massless quarks has been calculated in the mean-field 
approximation, which yields /XcCt ~ 0.66 [Hj. This analytic prediction has also been 
confirmed by numerical simulation of a gas of loops, the monomer-dimer-polymer 
ensemble, whose partition function is identical to that of lattice QCD with stag- 
gered fermions at /3 = ^21, and by the Glasgow reweighting method [13]. The 
numerical simulation of ^21 gives /Xc^ = 0.63(2) in the chiral limit and 0.690(15) 
at rriga = 0.1. The result from is /^cO = 0.687(15) at m^a = 0.1. At first sight, 
this excellent agreement between mean-field and Monte Carlo is very satisfying. 
On second thought, however, it is rather mysterious. One would expect the tran- 
sition to occur when (3/i) exceeds the free energy required to create an additional 
baryon, which is equal to its mass at T = 0. Now, every analytic prediction for 
the nucleon mass mjv gives a value (mjvct) very close to 3 so that {ficO') should be 
very close to 1. The values reported above are considerably smaller. Among the 
various explanations for this puzzle, we want to check the validity of the analytic 
approximations: is the nucleon mass really close to 3 in the strong coupling limit, 
or is it close to Sfic^ instead? 

In the next section (Sec. II), we review the analytic approximations which have 
been proposed for the hadron spectrum. Then (Sec. Ill) we present our simulation 
results and compare our numerical results to the analytic predictions in Table I. 
Conclusions follow. 
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2 Overview of analytic approaches 




In the strong coupling limit, the gauge field part of the action can be dropped 
because its coefficient is /3 = 0, and the gauge link variables become random SU{N) 
group elements. Then, gauge links give non-vanishing contributions only when a 
given link forms a part of an SU{N) singlet combination. Using this property, 
one can construct either an effective action and derive the hadron spectrum from 
it 12 El; or one can count the graphically relevant diagrams and obtain from them 
the hadron spectrum 

In the authors derive the effective action for SU (3) in terms of meson and 
baryon degrees of freedom by integrating the gauge link variables out. Then, via 
systematic large-ci (the spacetime dimension) expansion of this effective action, the 
chiral condensate and the correlation functions of the composite mesonic and bary- 
onic fields are obtained. The positions of poles in the zero momentum projected 
correlators give meson (Mm) and baryon (Mb) masses which can be summarized 
as follows: 

cosh M^a = d{T - 1) + 2A; + 1, (2) 
sinh Msa = ^A^(2rf)^ (3) 

where A = m + + 1, fn = ma/\/2d, ma is the dimensionless quark mass, 
and k = 0, 1,2,3 gives access to four meson channels, to which are assigned the 
physical meaning of the vr, p, oi, ao//o respectively. For small ma <^ V2d, 

ii^ip) ^ 2.12-0.75ma (4) 

{M^af ^ 5Mma (5) 

Mpa ^ 1.76 + ma (6) 

Ma,a ^ 2.29 + 0.58ma (7) 

M^^/f^a ^ 2.63 + 0.41ma (8) 

Mea ^ 3.12 + 1.06ma (9) 

Leading 1/d corrections were obtained in [S], yielding for d = 4, N = 3: 

(^tp) ^ 1.99 -0.56ma (10) 
(M^a)^ ^ 4.60ma (11) 
Mea ^ 2.93 + 1.99ma (12) 

and the masses of the other 3 mesons are unchanged. Leading corrections in 
were also obtained in Ref. [Hj. They always come in the combination l/{g'^N). 



3 



Perhaps this motivated the authors of |5 to try a different approach. They 
developed a graphical method for summing relevant, tree-like diagrams in the 
large-A^ limit {g'^N oo) at fixed d. They showed 

1 — dm^ — {2d — l)ma 

iV^^^^ = d- + (may 

cosh Mmc = ^ — + 2k, (14) 

zrrir 

sinh M^a = 2^^^ ^^^^ 



where = ma + y {may + 2d ~ 1 and k = 0,1, 2, 3. Again, for small ma <C 
^/2d — 1, one obtains 

{M^af ^ 4.54ma (16) 

Mpa ^ 1.76 + 0.81ma (17) 

Ma,a ^ 2.63 + 0.32ma (18) 

Mao/foa ^ 2.29 + 0.46ma (19) 

Mea ^ 2.91 + 1.14ma (20) 

The three sets of predictions are collected in Table 1 below, where they are 
denoted respectively "M.F."(mean field) [HIS], "l/d-corr." ^ and "d = 4, large- 
iV" 13] . The last two columns show our Monte Carlo results, presented in the next 
section, for the quenched and the full QCD theories. 

3 Numerical results 

We have simulated the partition function Z = J VU e^^ , with the staggered action 
S: 

S = Y1[^Y1 Vf^{x){i^i.x)U^i{x)ij{x + fi)- i){x + fi)Ui,\x)i){x)} + mail)il){x)] (21) 

X ^ ^1 

with rip{x) = (^_)^i+- +^M-i^ for 3 values of the quark mass {ma = 0.1,0.05 and 
0.025), on an 8^ x 16 lattice. In each case, 500 decorrelated configurations on 
which we measured the quark condensate and the hadron correlators have been ac- 
cumulated. Since, with staggered fermions, a given hadron correlator also contains 
contributions from its parity partner, we perform joint fits of the two correlators 
of parity partners 



C^{t) = A(e-'"^* + e-"^^(^-*)) + (-l)*fi(e-™^* + e~'"^(^~*)) 

C2{t) = A'(e-™2* + e-'"^(^-*)) + (-l)*5'(e-"^* + e-'"^(^"*)). (22) 
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Table 1: Comparison of hadron properties. 





rriqa 


M.F. 


1 / d-COTT. 


d = 4: 


MC results 


MC results 










large- H] 


quenched 


Nf = A 


average 


0.1 


— 


— 


— 


0.00001(4) 


0.00603(2) 


plaquette 


0.05 


— 


— 


— 


same 


0.00656(3) 




0.025 


— 


— 


— 


same 


0.00682(3) 




0.1 


2.0463 


1.9325 


1.9281 


1.9295(5) 


1.9172(2) 




0.05 


2.0838 


1.9606 


1.9562 


1.9574(7) 


1.9413(3) 




0.025 


2.1026 


1.9747 


1.9703 


1.9710(11) 


1.9520(6) 




0.1 


0.7521 


0.6780 


0.6735 


0.6776(1) 


0.6780(1) 




0.05 


0.5318 


0.4794 


0.4762 


0.4778(1) 


0.4784(1) 




0.025 


0.3761 


0.3390 


0.3367 


0.3374(1) 


0.3379(1) 




0.1 


1.863 


1.863 


1.841 


1.847(2) 


1.831(2) 


nip 


0.05 


1.813 


1.813 


1.801 


1.802(4) 


1.774(4) 




0.025 


1.788 


1.788 


1.780 


1.768(6) 


1.779(7) 




0.1 


2.350 


2.350 


2.336 


2.337(6) 


2.276(6) 




0.05 


2.321 


2.321 


2.313 


2.317(9) 


2.274(10) 




0.025 


2.307 


2.307 


2.302 


2.325(18) 


2.300(19) 




0.1 


2.675 


2.675 


2.662 


2.637(16) 


2.511(23) 




0.05 


2.654 


2.654 


2.646 


2.588(27) 


2.551(27) 




0.025 


2.644 


2.644 


2.638 


2.598(51) 


2.636(47) 




0.1 


3.225 


3.129 


3.024 


2.961(3) 


2.931(3) 




0.05 


3.172 


3.030 


2.967 


2.890(4) 


2.863(5) 




0.025 


3.146 


2.980 


2.939 


2.883(8) 


2.831(10) 



For example, the vr channel correlator is fitted simultaneously with the scalar 
channel correlator, and the p channel correlator with the pseudovector channel 
correlator by use of the 6-parameter fitting form of Eq. ()22|) . 

A few words on our simulation algorithm may be useful. The 4-flavor the- 
ory can be simulated with standard Hybrid Monte Carlo. Nevertheless, we used 
RHMC pn], with two sets of pseudo-fermion fields each yielding ^det(M) after 
integration. This decomposition allowed us to use a larger stepsize ^1]. More 
importantly, the stopping criterion of our solver was set to the very loose value 
1.0 X 10^^, while maintaining an acceptance ~ 65%. This great advantage of 
RHMC, recognized in ^Hl; provides an important gain in efficiency. 

Our spatial volume is reasonably large compared to the pion mass {irtT^L = 2.6- 
5.3), and is very large compared to the nucleon mass {L ~ 5.5 fm). Thus we expect 
negligible finite-size effects on all observables except perhaps the pion mass, and 
analytic predictions obtained on an infinite lattice El IH E] can be directly 
compared with our numerical results. Our m-,^/mp mass ratio ranges from 0.3703(5) 
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Figure 1: Hadron masses as a function of the quark mass m^a. The sohd hues 
show the predictions of Ref. j^. 



down to 0.1899(7), compared to the physical value 0.179. So we are studying the 
light-quark regime. Note, however, that mN/mp is about 1.6, close to the non- 
relativistic value 3/2, and hardly changes with the quark mass. So the spectrum 
of strong couphng QCD is only crudely similar to that of continuum QCD. 

Our results are collected in the rightmost two columns of Table I, for the 
quenched {Nj = 0) and the full {Nf = 4) theory. Ref. jlj counts all quark "tree 
graphs" which enclose zero area and is very similar to the quenched approximation. 
This motivated us to include quenched simulation results for comparison. Indeed, 
the tree graph resummation provides an excellent approximation to the quenched 
results^. Differences (of about 2%) are visible in the nucleon mass only. In turn, 

■^Compared to the quenched theory, mesonic tree graphs do not include the contribution of 
baryon loops, but these are suppressed by the heavy baryon mass. Baryonic tree graphs do not 
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Figure 2: The ratio of our Monte Carlo data over the analytic result of Ref. |3], as 
a function of the quark mass rriqa. Some channels are horizontally shifted to avoid 
overlapping data points. 



the quenched and full QCD results are very close to each other, with the most 
difference (at the 5% level) in the ao//o channel. 

The agreement between the full QCD results and the mean-field calculations [21 
0] is reasonable, and improves with 1/d corrections jH]. It is best in general with 
the large- approach of jlj. 

To illustrate the level of agreement with the predictions of |1] , we show in Fig. 1 
the masses we measured together with the analytic dependence of jl] as a function 
of the quark mass. In Fig. 2, we show the ratio of our measured values over those 
predicted by The agreement is better than within 1% for the chiral condensate 

include the contribution of meson loops, which is probably the reason why the nucleon mass 
prediction is slightly wrong. 
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2.0 



♦ mu_cfull model (Ref. 13) 
X monomer-dimer (Ref. 12) 
o m_baryon/3 (our data) 




CHEMICAL POTENTIAL 



Figure 3: mbaryonO./3 as a function of the quark mass m^a together with /icfl (our 
data is added to the Figure from Ref. [THj) 

{ipip) and the pion mass. The ratio for the rho mass varies non-monotonically, 
which is possibly caused by the crossing of the threshold for p — > vrvr decay. For 
the nucleon, the Monte Carlo data is within 5% of jl]. Our determination of the 
scalar mass (ao//o) is too poor to provide a real check. Interestingly, the lattice 
result agrees with the analytic result < nip < nia^ < maQ/f^, in contrast to the 
experimental ordering, m^r < rup < ma^/fQ < ma^. 

In Figure 3, we compare our nucleon mass with old Monte Carlo data for the 
critical, zero temperature quark chemical potential in infinite coupling QCD ^21 
IT^ . The error bar for 1/3 of our nucleon mass is smaller than the plot symbol. 
There is a clear, large difference between the critical quark chemical potential and 
1/3 of the nucleon mass. 
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4 Conclusions 



On the whole, our study justifies the mean- field approximation used in |21E1- It 
also provides a beautiful confirmation of the approach of |3], which goes beyond 
mean-field and becomes exact when g'^N ^ oo. It may seem strange in fact that 
this g'^N —>■ oo approximation is so good (at the percent level), since virtual quark 
loops are absent in this approximation in contrast to the numerical simulation. The 
explanation, we believe, lies in the small plaquette average value. The dynamical 
quarks shift the plaquette away from the value zero which it would take in the 
quenched case, but only slightly. This shift increases very mildly for lighter quarks. 
The lightest quarks which we considered, m^a = 0.025, give a pion to rho mass 
ratio near the physical one, but yield a plaquette value of 0.00682(3) only. This is 
much less than the naive estimate based on a 1/ (m^a) expansion ^21j , which would 
predict an effective shift A/5 ~ 0.8 in the gauge action, leading to a plaquette value 
(9(0.1), about 15 times larger than measured. Thus, the price to pay for creating 
virtual quark loops remains high, and their effect on the spectrum is small. For 
the same reason, the quenched approximation turns out to be closer to full QCD 
at strong coupling than at weak coupling: at strong coupling, the ordering effect 
of the fermion loops is "drowned" in the disorder of the gluons. 

As shown in Table I, the nucleon mass is just a bit smaller than 3. In particular, 
for rriqa = 0.1 which is the quark mass used in the finite density simulations 
of f^lE]; the nucleon mass is 2.931(3) and one would expect the T = transition 
induced by a quark chemical potential to occur at /ic = 2.931(3)/3 = 0.977(1). 
However, the measured value for fic |I2] is 0.690(15) (0.687(15) in [I2])- The 
discrepancy with the expected value is about 30%. In other words, the nucleon 
seems to have a free energy about 300 MeV less than its mass. Could this be caused 
by a strong attraction among nucleons ^j? The weakness of the real- world 
nuclear interaction results from near-cancellation between the attractive omega 
exchange and the repulsive sigma exchange fHl- This cancellation may not occur 
when the parameters of QCD are modified. Of course, other, less extraordinary 
causes may be at work: non-zero temperature corrections, or even algorithmic 
problems as suggested in This line of investigation should be pursued further. 
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